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Abstract 


The resporse of a vapor -gas bubble to a sinusoidal pressure pulsation Is 
analyzed. The esses of a spherical bubble in an infinite body of liquid and 
a truncated spherical bubble attached to a wall are examined. The Interior 
of the bubble is assumed to be a mixture of saturated vapor and noncondensable 
gas. The frequency response function is obtained and includes the effects of 
frequency, bubble size, and liquid subcooling. Experimental results of bub- 
bles in water and on metal surfaces are compared with the analysis. 


INTRODUCTION 

Two-phase flow dynamics has been a concern to the heat transfer engineer be- 
cause of the need to predict or to avoid boiling crisis. It has been shown 
(refs. 1,2) that premature burnout can result from flow instability in a 
"soft" system which is highly compressible. The response of void fraction to 
pressure perturbations is of great importance in analyzing two-phase flow sys- 
tem stability. 

The response of gas bubbles to pressure pulsations has been a subject of study 
in the acoustics field for many years (c.g., refs. 3,4). On the other hand, 
the behavior of vapor bubbles has only been studied for the case of bubble 
growth and collapse in a constant pressure field (refs. 5, 6, 7, 8). 

The purpose of this study is to examine the response of a vapor bubble in a 
saturated cr near saturated liquid to a pulsating pressure field. The effects 
of noncondensable gas, bubble shape and the presence of a metallic solid sur- 
face will be included in the general analysis. The resulting equations give 
the amplitude and phase lag of bubble size fluctuation in response to pressure 
fluctuation. The general equation can be reduced to special cases, such as 
those for a spherical bubble in a liquid bulk and a hemispherical bubble on a 
metal surface. 

Experiments were conducted to measure the bubble response to sinusoidal pres- 
sure pulsations in the above mentioned special cases. The experimental results 
are compared with the analysis. 

ANALYSIS 


Formulation 


Consider a spherical gas -vapor bubble in an infinite sea of motionless liquid. 
The geometry of the system is shown in figure l(a). The temperature of the 
liquid is assumed to be that for which a gas -vapor bubble of radius R Q with 
vapor partial pressure P yo and gas partial pressure Pg Q will exist in 
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equilibrium. If the pressure in the liquid begins to oscillate, the bubble 
size will oscillate because of the alternate superheating and subcooling of the 
vapor inside and the expansion and contraction of the gas. It will be assumed 
that during the oscillations the properties of the liquid remain constant, and 
that the interior of the bubble is homogeneous with temperature T b = T g = T y , 
pressure P b =» P y + P g and density * p y + p g . The liquid is assumed to be 
inviscid and incompressible. No heat sources are present. In such a case, the 
liquid momentum and continuity equations can be combined to yield the Rayleigh 
equation 


R^f ♦ 

dt^ 


3 dR 2 JL_ H. 
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The energy equation is 


(Sf + ^ • g™* hj 


r > R 


with boundary conditions 

Tj(R,t) - T fe 
Tj(r,t) bounded 

The equation of state for the gas is 

p g v % 

The vapor phase is assumed to satisfy the Claus ius-Clapeyron relation, 

h. 
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Combining equations (3) and (4) with T„ = Ty. = T b gives 
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dPg + dP y 
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The energy equation for the bubble interior is 
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where 


mC = m C + m C 

V v,v g g,v 

The terms on the left side represent respectively the energy increases in the 
bubble interior due to sensible heat, internal energy of vaporization, and 
compression of the gas and vapor. The term on the right side represents heat 
conduction from the liquid to the bubble interior. 


An expression for the rate of evaporation of the vapor must be obtained to 
complete the formulation. The ratio of masses of vapor ahd gas in the bubble 

is 
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Substitution of (5) and (8) into (6) yields 


( 7 ) 


( 8 ) 


— u„_m__P 

raC 

Py^b 


u« m p / n_ \ u_ m P 

+ ‘g v s (l + f e \ . f E v b 


v fg p g 


P V T b 


dT 


__b 

dt 


+ Ph(l - W) + 


u fg m v 1 dv 

v J dt 


dT, 

k l A b “^7* (**>*) ( 9 ) 


Solution 


The magnitudes of oscillations of all variables are assumed to be small enough 
that the system equations can be linearized by defining small dimensionless 
quantities x(t), 0 b (t), 0j(r,t), p(t), and c(t), each much smaller than 
unity, as follows: 


R = R 0 [l + x(t)], 

T z = T 0 [l + 0 j(r,t) ] , 

P b = P 0 tl + W + p(t)]. 
Dimensionless time and distance functions 


T b - Tot 1 + V*)! 
P« = P 0 Il + e(t)] \ 
W = 2a/B o P 0 
are defined as 


( 10 ) 


t = a z t/R 2 , y = r/R 0 


( 11 ) 


When the quantities defined in equations (10) and (ll) are substituted into 
(l), (2^, (5), and (9) and terms involving products of small quantities are 
neglected, the following equations result: 
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where 


B = P o R2/ 0l a| 



(15a) 


Equations (12) to ( 15 ) constitute a set of linear differential equations 
relating the dimensionless pressure variation €(t) to the dimensionless change 
in bubble radius x(t). It is shown in standard books on linear system theory 
that if such a system is excited by an oscillatory input, the response will be 
oscillatory. The ratio of the magnitudes of the output and input variables is 
given by the magnitude of the frequency response function, which is the ratio 
of the Laplace transforms of the output and input variables with zero initial 
conditions and with the Laplace variable S set equal to the imaginary dimen- 
sionless frequency 

jn = J^/ai) . 

Solution of equations ( 12 ) to (15) for the ratio £{x(t) ) /£{c (t) } = X(s)/e(s) 
is straightforward and the details will net be discussed here. The corre- 
sponding frequency response function is 


- I 

K 2 (l + W)jfl 

i + <\/jn + iqjn 

(16) 

Since X/E represents the ratio of the dimensionless radius change to the di- 
mensionless pressure variation, it can be called the dimensionless compress- 
ibility of a bubble. 
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Bubble on a Wall 


If the bubble rests against a solid wall (fig. 1(b)), heat conduction takes 
place through the liquid and also through the wall. The bubble is assumed to 
be a segment of a sphere, and equations (12) to (14) still hold true. The 
bubble energy equation (15) must be altered to include heat conducted through 
the solid wall. This equation becomes 



(17) 


The heat conduction equation for the solid wall must also be accounted for. 
In dimensionless form, 
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Boundary conditions are taken as 


0 w (o,t) 




(18a) 


When these changes are made the modified form of the frequency response func- 
tion is 
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When La^Q/R 0 u v is large, 
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EXPERIMENTAL APPARATUS AND PROCEDURE 


The experimental apparatus consists of a test cell, a pressure pulse gener- 
ating device, a camera, and a data recording device (fig. 2). 

The test cell is a cylindrical tank with an optically-flat front window. The 
cell is connected to the pulse generator .and the pressure transducer at the 
bottom. A vent is located at the top of the circular siuewall to eliminate 
vapor (or gas) pockets. A nylon net with fine mesh was used to trap the 
vapor bubbles for study of bubbles in the bulk. Metal blocks (about 2.5 cm 
inside), of copper and stainless steel were used for study of bubbles on 
metallic surfaces. The undersurfaces of these specimens were concaved 
slightly to trap bubbles. A 100-watt heater was situated under the surfaces 
to keep the water warm and to supply vapor bubbles. Additional heat was sup- 
plied by infrared lamps. The unexposed sides were insulated. 

The pulse generator was a four -way valve made to rotate. The four ports were 
connected to the laboratory low pressure air line, vacuum line, atmosphere 
and test cell (through a surge tank). The test cell was thus alternately ex- 
posed to positive and negative pressures, and a smooth sinusoidal pressure 
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wave was obtained. The frequency is adjustable through changing the rotating 

.pec 1 of the valve over a limited range of 1 to 4 cycles per second. 

The thermocouple readings of the bulk and the metal specimens and the trans- 
ducer readings were monitored by a digital voltmeter. Temperature cc^Jd be 
read to within 0.2° C accuracy. But due to the slight nonuniformity of tem- 
perature field in the cell, the accuracy of temperature data is probably no 
better than 1° C. Pressure could be read to within 1 cm of water column. 

For synchronization purposes, both the camera and the recorder were provided 
with 100 cps timing marks. A zero-time lamp was put in front of the c»-mera. 

When the lamp turned off a galvanometer reading would deflect on the oscillo- 

graph chart. These two events are taken as zero-time. 

During each run, while data was being taken, the heater was •'•urned off to en- 
sure clean sinusoidal pressure waves. Otherwise, the heater would generate 
a lurge number of bubbles, thus imposing random pressure pulses onto the p 2 *i- 
mary pressure wave. 


RESULTS AIJD DISCUSSION 
Presentation of Results 


Figures 3, 4, and 5 are comparisons between experimental and analytical re- 
sults for bubbles in the bull., on the copper surface, and on the stainless 
steel surface. Curves are plotted for the case R 0 = 10" ^ centimeter. The 
effect of radius is felt only at very high frequencies, however, and curves 
at other R^ values would be almost indistinguishable from those given in 
the frequency range shown. 


Controlling Mechanisms 

The analytical equations (16), (20) have incorporated the effects of inertia, 
heat transfer, and surface tension. Close examination of the relative magni- 
tudes of the various terms using properties of water reveals that the bubble 
response is inertia-controlled in the high frequency range (ft > lO) and heat 
transfer-controlled in the low frequency range (ft < 10 3 ). This is consistent 
with the bubble dynamic analysis for the constant pressure field (e.g. in 
nucleate boiling, cavitation, etc.) in which it is generally accepted that 
bubble growth rate is inertia-controlled only in the very early stage 
(t < 10“' J sec) but is heat transfer-controlled throughout most of the bubble 
life (t > 10“2 sec) . 

For the size of bubbles and the frequency range encountered in the experi- 
mental study, 0 is usually in the range of 10^ to 10 3 . In this range, heat 
transfer is the controlling mechanism. 


Effect of Frequency 

The qualitative trend of figures 3, 4, and 5 shows that when bubbles contain 
only vapor (Tj-, = 373° K for water under 1 atm), the amplitude response of the 
bubbles is strongly dependent upon the dimensionless frequency ft. These bub 
bles are highly compressible in the lower ft range and become far less com- 


7 


pressible at high dimensionless frequency. This behavior can be explained by 
consideration of the heat transfer mechanisms. The dimensionless frequency, 
by definition, is the ratio of the thermal diffusion time, P 2 /a to the pulsa- 
tion period. Therefore, when the dimensionless frequency is large, the time 
available during a compression or expansion period is insufficient for heat to 
diffuse through the liquid and the bubble behaves nearly adiabatically . In 
such a case the mass transfer due to latent heat exchange is small and the 
condensable vapor behaves similarly to the noncondensable gas. On the other 
hand, at lower frequencies case, a long period is available for the heat to 
diffuse. Thus, latent heat involved with condensation can be exchanged back 
and forth between the bubble and the liquid. The accompanying mass transfer 
gives rise to large bubble pulsations. When even a small amount of non- 
condensable gas is present in the bubble, the dependence of compressibility 
on frequency is greatly diminished. 


Effect, of Noncondensable Gas 


Th'* effect of noncondensable gas is shown by the lines marked with bulk tem- 
perature. It was observed that noncol lapsable bubbles could be obtained and 
maintained indefinitely in subcooled water. This indicates the presence of 
noncondensable gas in the bubble in addition to the vapor. Assuming that the 
partial pressure of the vapor can be expressed by the vapor pressure along 
the saturation curve, the balance of the partial pressure must be contributed 
by gas dissolved in the water which was released into the bubble. Thus, a 
bubble in the subcooled liquid contains a mixture of vapor and gas, with the 
composition of gas increasing with increasing subcooling. The analytical 
lines showed that a trace amount of gas (or a very small subcooling), could 
mean an order of magnitude drop in amplitude response. The experimental data 
for near saturation condition lay in between the saturation curve and 1° C 
subcooling curve. While the 2 ° C subcooled data lay Just below the corre- 
sponding 271° K line. This result tends to substantiate the analysis. It 
also chows how much one could be in error if a trace amount of noncondensable 
gas is discounted as negligible. It should be emphasized that in the present 
study, the bubble interior is assumed to be homogeneous. Thus, retardation 
of diffusion of vapor caused by the noncondensable gas is neglected. This is 
Justified on the grounds that mass diffusion due to the change of phase for 
a pulsating bubble takes place in alternate directions, and thus it would not 
cause an accumulation of inert gas to form a diffusion barrier as observed in 
the film condensation studies (ref. 9 ). 


Effect of Surface 


The presence of a solid surface provides another avenue for heat to be trans- 
ferred in and out of the bubble. In this analysis, it is assumed that the 
liquid film between the bubble and the solid surface is so thin that it does 
not offer any resistance to heat transfer. Such an assumption might be vul- 
nerable for the case of high frequency. For the low frequency range, the 
experimental results appear to agree with the analysis reasonably well. It 
should be noted that, in this study, thick blocks of metals were used. There- 
fore, only the effect of conductivity and thermal diffusivity can be studied. 
To this effect, the results show that bubbles on copper are more compressible 
than stainless steel which is less conducting. Equation ( 19 ) indicates (with- 
out experimental verification) that the thickness of the metal wall could be 
a factor too. Thus, the void response in a tube to a pressure pulse could 
very well depend upon the tube thickness as well as the tube material. 
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The preaenee or metal surface effects the dimensionless compressibility more 
strongly when the liquid temperature approaches saturation temperature (l,e., 
near al-1 -vapor condition). The effect of a solid surface diminishes as gas 
content to increased with increasing subcooling. 


Effect of Bubble Shape 

The bubble shape does not appear to greatly alter the compressibility . Ap- 
parently due to the fact that the ratio v/A^ does not vary over a very great 
rangu (o.g. , for the same volume, the surface areas of spherical bubble and a 
hemispherical bubble differ only by 20 percent). Thus, for estimation pur- 
pose, a spherical bubble could be assumed without introducing a large error. 


C ONCLUSION 

The dynamic response of a vapor bubble to a sinusoidal pressure pulsation was 
studied both analytically and experimental.] y, for the cases of bubbles in the 
liquid bulb and on metal surfaces. The studies showed that bubble response is 
strongly dependent on the frequency and gas content. It is also dependent on 
the thermal properties of the metal but is only weakly dependent on the bubble 
shape. The analysis also showed that the bubble dynamics are inertia con- 
trolled in the high frequency range but thermal -controlled in the low fre- 
quency range. 


NOMENCLATURE 


A 

B 

C 

E 

h fg 


d 


1 


L 

M 

m 

P 

P 

R 

r 


area 

a dimensionless parameter defined in equation (15a) 
specific heat 

Ap/pQ 

latent heat of evaporation 
thermal conductivity 

dimensionless parameter defined in equation (15a) 
dimensionless parameter defined in equation (l5a) 

thickness of solid wall 
molecular weight 
mass 

pressure 

small perturbation in pressure, dimensionless, equation (10) 
bubble radius 

coordinate as shown in figure 1 
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T 

t 

u 

? 

V 

v fg 

W 

X 


y 

z 


a 

e 

6 

a 

P 

i 

T 

ft 


temperature 

time 

internal energy 
velocity 

volume 

change of specific volume due to evaporation 
Weber number, defined in equation (10) 

ar/r 0 

small perturbation in radius, dimensionless, equation (10) 
dimensionless radius coordinate, equation ( 11 ) 
compressibility coefficient 
coordinate, figure l(b) 
thermal diffusivity 

small perturbation in pressure, dimensionless, equation (lO) 
small perturbation in temperature, dimensionless, equation (lO) 
surface tension 
density 

dimensionless coordinate, z/R 0 

dimensionless time, equation (ll) 
dimensionless frequency, 2 ti:/t 


Subscripts: 
b bubble 

fg change from liquid to vapor 

g gas 

l liquid 

o unperturbed value 

v vapor, also constant volume 

w wall 


infinity 
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